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Abstract 

We study the structure of 2D electronic states in a strong magnetic field 

in the presence of a large number of resonant scatterers. For an electron in 

the lowest Landau level, we derive the exact density of states by mapping 

the problem onto a zero-dimensional field-theoretical model. We demonstrate 

that the interplay between resonant and non-resonant scattering leads to a 

non-analytic energy dependence of the electron Green function. In particular, 

for strong resonant scattering the density of states develops a gap in a finite 

energy interval. The shape of the Landau level is shown to be very sensitive 

to the distribution of resonant scatterers. 
Pacs numbers: 73.20.Dx, 73.40.Hm 
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During recent years there has been a growing interest in the role of multiple resonant 
scattering in transport. Most of the studies have been related to the passage of light through 
a disordered medium. In particular, it was shown in a recent experiment |]J and subsequent 
works [f2] , that multiple scattering near resonances leads to a renormalization of the diffusion 
coefficient up to an order of magnitude. 

It is natural to expect that resonant scattering would also affect quite strongly the 
properties of electrons in disordered systems. The effective trapping of the electron in 
resonant states is expected to suppress diffusion, just as in optics fTU^. This would in turn 
be evident in the single-particle density of states (DOS) and localization properties. 

In this paper we study the electronic states of a 2D system in a strong magnetic field in 
the presence of a large number of resonant scatterers. This choice is motivated in part as 
experimental structures with such geometry became recently available thanks to remarkable 
advances in the fabrication of arrays of ultra-small self-assembled quantum dots jj. With 
typical sizes of less than 20 nm and very narrow variations of less than 10%, an array of 
such dots with density 10 10 — 10 11 cm" 2 can be produced at some preset distance from a 
plane of a high mobility electron gas. Q As the Fermi energy in the plane approaches the 
levels of dots, the virtual transitions between dots and the plane result in multiple resonant 
scattering. Such scattering which, in principle, extends through the entire system, strongly 
affects the DOS of a 2D electron. 

The DOS of 2D disordered electronic systems in a quantizing magnetic field has been 
extensively studied for the last two decades |5|-[T5| . The macroscopic degeneracy of the Lan- 
dau levels (LL) makes impossible a perturbative treatment of even weak disorder and calls 
for non-perturbative approaches. For high LL, Ando's self-consistent Born approximation 
H was shown to be asymptotically exact for short-range disorder [|IT| , |I3|| , while in the case 



of long-range disorder the DOS can be obtained within the eikonal approximation |13[ . For 
low LL and uncorrelated disorder, the problem contains no small parameter and neither of 
those approximations apply. Nevertheless, Wegner was able to obtain the exact DOS in a 
white-noise potential for the lowest LL, by mapping the problem onto that of the 0D com- 



plex 4 -model 0. This remarkable result was extended to non-Gaussian random potentials 



by Brezin et al. PH |, and recently to multilayer systems 

The "regular" disorder broadens the LL into a band of width T. At the same time, the 
resonant scattering leads to a sharp energy dependence of the DOS near the resonance. The 
scattering is enhanced close to the LL center and is suppressed in the tails. Therefore, the 
efficiency of resonant scattering is characterized by the ratio 7/r, where 7 is the width of 
the energy spread of resonant states. 

The interplay of the resonant and non-resonant scattering leads to a rather complex 
energy dependence of the DOS. Nevertheless, for the lowest LL the problem can be solved 
exactly [see Eqs. (||) and (J9|) below]. We exploit the hidden supersymmetry of the lowest LL 
PJlO[| in order to map the averaged Green function onto a version of OD field theory. The 
DOS appears to be non-analytic as a function of energy; in particular, it develops a gap as 
resonant scattering becomes strong. 

The model. — Consider a 2D electron gas separated by a tunneling barrier from a system 
of localized states (LS). In addition to LS, a Gaussian random potential V(r) with correlator 
(V(r)V(r')) = w8(r — r') is present in the plane. We assume that energies of LS are close 
to the lowest LL and adopt the tunneling Hamiltonian 

where e M , and c M are the eigenenergy, creation and annihilation operators of the eigenstate 
of the Hamiltonian H + V(r) (H describes a free electron in magnetic field), e^, cj and 
Cj are those of the ith LS, and t^ is a tunneling matrix element. The latter is defined as 
= f drdzip*(r, z)Vi(r, z)ipi(r, z) ~ V^( r i, z d I drdzV^r, z)^(r, z), where V^r, z) is the LS 
potential and ipi(r,z) is its wave function. In the direction normal to the plane, the wave 
function ^(r, z) decays as e~ KZ , while in the plane it behaves as an eigenfunction ^*(r) of 
the Hamiltonian Hq + V(r). For high enough tunneling barrier, the dependence of k on /1 
can be neglected |TB[, so that t^ ~ %fj*^{vi)ti. 



A formal expression for the Green function of a 2D electron with energy u, G^ u (oj 



(fi\(uj — H) can be derived by integrating out the LS degrees of freedom. It has the 

form G(uj) = (u — e — S) -1 , where e is a diagonal matrix with elements e M , and self-energy 



matrix, 



E 



(2) 
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comes from scattering of the electron by LS. In such a form, however, the Green function is 
hard to analyze. Instead, it is convenient to work with an effective in-plane Hamiltonian, 
H ef[ , for the electron with energy uj. Recasting E Mi/ in coordinate representation, we obtain 
H e s(id) = Hq + V(r) + U(u, r), where the last term, 



describes the resonant scattering of electron by the LS. The potential (||) resembles that 
of point-like scatterers. The crucial difference, however, is that here scattering strength 
depends on the proximity of the electron energy to the LS levels. It is important to notice 
that U(u, r) changes from repulsive to attractive as the electron energy passes through the 
resonance. Since positions of LS are random with uniform density n LS , the distribution 
function of U is Poissonian. 

In the following, we assume that the tunneling barrier is high enough, and neglect the 
difference in ti for different LS, setting ti — t in the rest of the paper. Strong magnetic field 
implies that scattering keeps electron in the lowest LL. While for the white-noise potential 
this condition is standard, it is more restrictive for the resonant scattering. It should be 
noted, however, that the latter is effectively reduced by the energy spread of LS. 

The calculation of the DOS, g(id) = —7i^ 1 lmG(r,r), requires averaging of the Green 
function, G(r, r) = (r|(aL|_ — H c g)~ l \r) (with ui + = id + iO) over both random potentials V(r) 
and U(uj). Below, we derive this DOS exactly by using the approach of Ref. ||10|| . 

Derivation of DOS. — The Green function is presented as a bosonic functional integral 
G(r, r) = — iZ^ 1 j VipV%[)e lS %lj(v)%l)(r) with the action S[?p,ip} = J drip(r)[id + — H e g(uj)]i/;(r). 
After writing the normalization factor as a fermionic integral Z~ x = J DxDx e%S with the 
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u(u,t)=y: 



Id - 6i 



S(*i - r) 



(3) 



same action S[x, x]> both ip and x are projected on the lowest LL subspace as (u}—H )ip = uip 
(measuring all energies from the lowest LL). In the symmetric gauge, this projection is 
achieved with ip = (27r/ 2 )~ 1//2 ''e~^ ''> ' u{z) and x — (27r/ 2 ) _1 / 2 e~l 2 l 2//4i2 t>(z), where the bosonic 
field u(z) and the fermionic field v(z) are analytic functions of the complex coordinate 
z = x + iy (I is the magnetic length). The Green function then takes the form G(r, r) = 
— i(27r/ 2 ) -1 e - ! 2 ! 2 / 2i2 (u(z)u(z*)) , where (• • •) denotes a functional integral over u(z) and v (z) 
with the action 

S = J ^- {zl2/2l2 (uu + vv)[u + -V- U(u)). (4) 

As a next step, one introduces Grassman coordinates 9 and 9*, normalized as 
/ d 2 zd 2 9e~^ 2 ~ ee * = 1, and defines analytic "superfields" &(z,9) = u(z) + 9v(z)/y/2l and 
$(z*,9*) = u(z*) + 9*v(z*)/\^2l, taking values in the "superspace" £ = (z, 9). Using 
(u) = (v ) = and (uu) = (vv), the Green function can be presented as 



.e 



-er/a 
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G = J V$V$e tb $(£)$(C), (5) 

where ££* = \z\ 2 + 99* and S[Q>, $] is obtained from (f|) by substituting uu + vv = 

27r/ 2 /rf 2 ^e- er / 2p l>(r)$(0- 

We now perform the ensemble averaging over V and U. The Gaussian averaging of 

expiJVQd 2 z, where Q = J d 2 9e~tt* / 2l2 gives exp [-(w/2) J Q 2 d 2 z], while the 

averaging of expz J UQd 2 z with a Poissonian distribution function yields [17| 

^^VA (6) 



exp< -n l& 



1 — ( exp 



u — e, 

where (• • -) e denotes energy averaging. As a result, one obtains the following effective action 

r 2 r d 2 z < - n 2 



$] = tu; + J d 2 £a - — J |A | d 2 fl a 



/ 1 1 - ( ex P 



(7) 



2tt/ 

where a(f,f*) = e~^*^ Here T = (w/2^1 2 ) 1 ' 2 is Wegner's width of lowest LL 
(in the absence of resonant scattering), v = 2nl 2 n LS is the "filling factor" of LS, and we 
denoted A = 5 2 /(cu — e), where 5 = £/(27r/ 2 ) 1//2 characterizes the tunneling, 



The action (^) possesses a supersymmetry, characteristic for the lowest LL 0|lOj- Being 
evident for the first term, this symmetry between z and 9 can be made explicit for the sec- 
ond and third terms also by making use of the identity |10] n (2,7rl 2 J d 2 9e~ ee */ 2l2 <b&J = 
2irl 2 Jd 2 9e- n99 */ 212 which allows one to replace any functional of the form 

/ d 2 zf [2-kI 2 J d 2 9a) with 2tt1 2 J d 2 £h(a), where dh(x)/dx = f(x)/x. As a result one ob- 
tains a manifestly supersymmetric action S = J d 2 ^A(a), where 

T 2 a 2 



iA(a) = ioo + a 



-v 
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The supersymmetry leads, in turn, to the exact cancelation of contributions from z and 9 
spatial integrals into each diagram, so that the entire perturbation series can be generated 
in the OD field theory with the same action The Green function is then given by the 

ratio of two ordinary integrals, G(uS) = — i(2irl 2 )~ 1 ZQ J d 2 (j)e lA (jxj)*, where Z = J d^cpe^ 
with A(4>4>*) from (§). From this Green function, the DOS is obtained as 



g{u) = In r~dae^\ (9) 



2n 2 l 2 du + Jo 
where the derivative applies only to the first term of (H). 

Examples. — The energy averaging in @ can be performed analytically for an arbitrary 
distribution of LS levels, / 7 (e — e), where e is average energy and 7 is the width. The result 
reads 

T 2 a 2 



iA(a)= ioo + a 
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•v I" — / 7 (x)e i( "- e > [l - J (26y/xa)] , (10) 
J *c 
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where / 7 (x) is Fourier transform of / 7 (e) and Jo is the Bessel function. Numerical results 
for DOS with Gaussian distribution, / 7 (x) = e _7a;2//2 , are presented in Fig. 1. 

Consider first the case of a strong in-plane disorder, T/5 ^> 1. For a not very small 7, so 
that 5 2 /^Y <C 1, the Bessel function in ( ]Tt]| ) can be expanded to first order, yielding G{uS) = 
GwiyJ ~ S), where Gwiyj) is Wegner's Green function (that is with v = 0) and £(u;) = 



—iv5 2 Jq 00 cfee -7 x / 2 + l ( w - e ) x i s the first-order self-energy due to the resonant scattering. If 
the resonant level is close to the LL center, uj ~ e <C T, the first-order correction to the DOS 
reads 



6g(u) 7r — 2 v5" {'jj - c 



exp 



2 7 2 



2" 

(11) 



g W (0) y/2 7F 

where gw{w) is Wegner's DOS. 

Resonant scattering in this case manifests itself as a minimum of width 7 on top of the 
wider peak of width T. The evolution of the DOS with increasing S/'y is shown in Fig. 1(a). 
The effect is strongest for 5/7 ^> 1, however splitting remains considerable even for 7/5 ~ 1. 
For 5/7 1 the DOS is basically unaffected by resonant scattering and reduces to Wegner's 
form gw(u). 

With increasing scattering S/T, the shape of the DOS undergoes drastical transformation 
[see Fig. 1(b)]. For a strong scattering, the DOS develops a gap in the energy interval 
00(00 — e) < 0. The existence of the gap can be traced directly to Eq. (§) (with vanishing 
7/5 and T/8). In this energy interval the integration path in the a-integral in (^) can be 
rotated by e - i7rs s n ^- e ")/ 2 , resulting in a purely real iA. The origin of the gap is the following. 
If the "regular" disorder is weak (small T), the LL broadening comes from the resonant 
scattering alone. Then the scattering potential @ appears to be attractive for u < e, 
pulling the electronic states from the LL center to the left, while for uj > e the potential is 
repulsive, pushing the states to the right. Note that for a low density of scatterers, v < 1, 
a fraction 1 — v of states in the plane remains unaffected. Such "condensation of states" 
was known also for the case of repulsive point-like scatterers with a constant scattering 
strength [^,0.|lO|.[lTJ . In fact, the analogy extends also to the intricate structure of the DOS 



away from the gap. In particular, the smaller peaks correspond to singularities in g{uS) at 
integer values of oj{uj — e)/<5 2 [fT0[| . The behavior of g(uS) near the gap edges is different for 
uj — > and uj —>■ e: one can show that in the former case the DOS exhibits a discontinuity, 
g{uj) oc (1 — v)5(uj) + const/ 1^1", while near the resonance it vanishes as {uj — e) 1 ^" '. With 
increasing 7, the gap and the smaller peaks are washed out; however the peak at uj = 



persists throughout [see Fig. 1(c)]. 

In conclusion, although our derivation was restricted to the lowest LL, we believe that 
our results are more general and valid for higher LL also. Indeed, the gap in the DOS for 
small disorder is apparently a result of the LL degeneracy. Therefore, the above argument, 
related to the change in the sign of the potential (|3|), should hold for arbitrary LL. Note 
that the "condensation of states" also occurs for all LL numbers [fLT| . Thus, we expect that 



the gap in the DOS will persist, although the precise behavior of g(u) near the gap edges 
could be different. Concerning the sharp minimum in the DOS in the absence of the LL 
degeneracy [see Fig. 1(a)], it seems that this is a rather general feature. In fact, in the 
absence of magnetic field, analogous behavior has been known in the 3D case for identical 
scatterers 

A possible experimental realization of the multiple resonant scattering could be a system 
of self-assembled quantum dots separated from a 2D electron gas by a tunable tunneling 
barrier Due to the narrow distribution of dots' sizes, the spread in their energy levels, 7, 
does not exceed 10 meV ||. For considerable effect of the resonant scattering, one must have 
5 2 /^T ~ 1. For a typical LL width T ~ 1 meV, this condition implies that the parameter S 
should be about several meV, which would be reasonable to achieve. Moreover, for S/T ^ 1, 
an even weaker condition, the tunneling would be relatively strong and the effect of the 
resonant scattering would be significant. It was observed in Ref. |4j] that the mobility of the 
2D gas (at zero magnetic field) dropped by two orders of magnitude when the thickness of 
tunneling barrier between the dots and the plane was reduced. Although, we cannot give 
quantitative estimate for the zero-field case, this is certainly in qualitative agreement with 
our results. We hope that our results would further motivate experiments in magnetic fields. 

Finally, we have disregarded the possible charging effects and assumed that the tran- 
sitions occur between the plane and unoccupied dots. Certainly, as the Fermi energy ap- 
proaches e some of the dots will become singly occupied. Once occupied, such dots would 
have much higher energies and would not participate in the resonant scattering, reducing 
the effective density n LS , apart from producing (uniform) Coulomb shifts in the energies of 

8 



unoccupied dots. 

This work was supported in part by the US Department of Energy Grant No. DE-FG02- 
91ER45334. 
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FIGURES 

FIG. 1. (a) DOS [in units of g x = (2nl 2 )- 1 T- 1 ] for strong in-plane disorder, 8/Y = 0.3, with 
e = and v = 1.5, is shown for different 7/(5 = 0.1 (solid line), 0.5 (dotted), 1.0 (dashed), 2.0 
(long-dashed), and 10.0 (dot-dashed), (b) DOS [in units of 52 = {2tt1 2 )~ 1 S~ 1 ] for strong tunneling, 
<J/ 7 = 10.0, with e = 5 and v = 0.8, is shown for r/<5 = 0.1, (solid line), 0.2 (dotted), 0.3 (dashed), 
0.5 (long-dashed), and 1.0 (dot-dashed), (c) The DOS for weak in-plane disorder, F/S = 0.1, 
with e = 5 and v = 0.8, is shown for 7/T = 1.0, (solid line), 3.0 (dotted), 5.0 (dashed), and 10.0 
(long-dashed) . 
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